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instrumented  with  thermocouples,  heat  flow  meters,  and  calibrated  gutters.  Measurements  were  recorded  for  the 
winters  of  197 1  -72  and  1972-73.  The  degree  of  icing  and  the  chronological  changes  in  the  snow  cover  were  recorded 
on  35-mm  Kodachrome  slides.  It  was  found  that  eave  icing  is  a  sensitive  function  of  the  slope,  roof  covering  composi¬ 
tion,  and  solar  radiation.  The  effects  of  wind  were  not  investigated;  the  data  were  screened  to  remove  all  informa¬ 
tion  corresponding  to  windspeeds  over  8  km/h.  In  order  of  increasing  tendency  to  form  ice  dams  on  the  eaves,  the  _ 
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SUMMARY 


The  response  of  a  roof  to  its  snow  cover  and  the  environment  is  a  complex  phenomenon  de¬ 
pending  sensitively  on  the  design  and  construction  of  the  roof.  This  study  investigated  the  roles 
of  the  slope  and  composition  of  the  roof  covering. 

When  there  is  no  snow  cover  before  snow  begins  to  fall,  the  slope  of  the  roof  dominates  its 
icing  behavior  by  controlling  the  interaction  with  solar  radiation  and  attic  heat,  the  rate  of  fall  of 
snow  sliding  from  the  roof,  and  the  rate  of  flow  of  meltwater  from  the  roof.  In  general,  as  the 
slope  increases,  the  attic  air,  and  the  middle  section,  eaves,  and  lower  section  of  the  roof  increase 
in  temperature  while  the  temperature  of  the  ridge  decreases. 

Interfering  with  the  action  of  the  slope  on  a  roofs  snow  cover  are  the  smoothness  of  the  roof 
covering  and  its  thermal  characteristics.  Shingles  with  large  butt  ends  impede  the  fall  of  snow  from 
the  roof.  If  its  solar  absorptivity  is  low,  the  roofs  surface  does  not  interact  strongly  with  solar 
radiation  even  though  optimally  exposed  by  a  high  slope.  If  the  roofs  thermal  conductivity  and 
diffusivity  are  low.  it  may  be  able  to  absorb  large  amounts  of  thermal  energy  without  a  change  in 
temperature  significant  enough  to  melt  its  snow  cover,  thereby  reducing  the  effects  of  solar  radia¬ 
tion,  attic  heat,  and  heat  trapped  under  the  eaves. 

Once  any  ice  has  formed  on  the  eaves,  it  acts  as  a  barrier  to  the  flow  of  meltwater  and  sliding 
snow  while  also  acting  to  maintain  the  eaves  at  a  cooler  temperature  than  the  middle  section  of  the 
roof,  producing  and  maintaining  the  temperature  profile  that  leads  to  more  icing.  Under  these 
circumstances,  the  roles  of  the  thermal  and  optical  characteristics  of  the  roof  dominate  the  icing 
behavior. 

To  reduce  or  eliminate  icing,  a  roof  may  be  designed  to  have  a  high  enough  slope  to  eliminate 
as  much  of  the  snow  cover  as  possible  through  the  action  of  gravity.  This  must  not  enable  the 
roof  to  interact  with  solar  radiation  and  attic  heat  enough  to  cause  a  temperature  rise  of  the 
middle  section  sufficient  to  generate  meltwater  that  may  become  ice  at  the  eaves.  The  thermal 
properties  of  the  roof  covering  must  be  chosen  to  compensate  for  higher  slope.  An  alternative 
approach  is  to  design  the  roof  with  a  low  slope  and  rough  surface  so  that  its  exposure  to  solar 
heating  is  reduced  and  a  sufficiently  deep  snow  cover  is  retained  to  insulate  its  surface  from  the 
exposure  to  solar  radiation  that  cannot  be  avoided.  In  the  latter  case,  it  is  important  to  insulate 
the  attic  floor  and  the  underside  of  the  eaves  while  constructing  the  roof  of  materials  of  low 
thermal  conductivity. 

In  order  to  promote  improvements  in  roof  design,  it  would  be  wise  to  direct  future  research 
into  better  measurements  of  the  thermal  and  optical  properties  of  roof  materials  and  to  study  the 
detailed  interactions  between  snow  cover,  roof  characteristics,  and  environmental  variables. 

Efforts  are  currently  being  made  in  these  directions.'  ‘ 


vi 


ROOF  RESPONSE  TO  ICING  CONDITIONS 


J.W.  Lane,  S.J.  Marshall  and  R.H.  Munis 


INTRODUCTION 

Roofs  for  use  in  cold  regions  should  be  designed 
to  minimize  icing  and  enhance  the  thermal  integrity 
of  the  total  building.  The  former  problem  has  long 
been  considered  important,  and  the  recent  empha¬ 
sis  on  energy  conservation  has  focused  attention  on 
the  latter.'  *  The  purpose  of  this  work  is  to  des¬ 
cribe  the  thermal  response  of  a  roof  (as  a  function  of 
its  design)  to  the  conditions  that  promote  icing. 

The  variables  with  which  this  work  will  deal  can 
be  divided  into  three  different  types:  1)  variables 
that  specify  the  design  of  the  roof,  “roof  character¬ 
istics”;  2)  variables  that  specify  the  thermal  response 
of  the  roof,  “roof  responses”;  and  3)  variables  that 
specify  the  icing  conditions,  “meteorological  vari¬ 
ables.”  Table  I  lists  these  variables  and  defines  them. 
As  others  could  have  been  chosen,  this  list  serves  to 
define  the  scope  of  this  discussion. 


heat  flow  profile  of  /th  roof 
snow  depth  of  /th  roof 
%  of  area  covered  at  depth  hj 
volume  of  meltwater  runoff  for  /th  roof 
degree  of  icing  for  /th  roof 
1  means  solid  layer  of  ice  under  snow 
0.75  means  solid  layer  of  ice  beginning  to 
form 

O.SO  means  numerous  icicles  extending 
from  snow  cover  or  roof  covering 
0.25  means  few  small  icicles  hanging  from 
eave 

0  means  no  visible  sign  of  ice 

Meteorological  variables 

vertical  eppley  solar  radiation 
S  snowfall 

t  time  of  day 

7X  dir  temperature 
pj  density  of  snow 

density  of  water 
V  windspeed 

EXPERIMENTAL  PROCEDURE 


"j 

Tj 

t'i 


Table  1.  Definition  of  variables. 

Roof  characteristics 

A  thermal  diffusivity  of  plywood  sheathing 

A*  thermal  diffusivity  of  shingles 

c  composition  -  asphalt,  aluminum,  cedar 
i  section  of  roof  (E  =  eave,  L  -  low,  M  -  middle, 

R  =  ridge,  U  -  underdeck,  F  =  attic  floor) 

/  (a,  c')  =  1  to  4  (1  =  low-slope,  asphalt  roof;  2  =  low> 

slope,  aluminum  roof;  3  =  low-slope,  cedar  roof; 

4  =  high-slope,  asphalt  roof;  5  •  high-slope,  alumi¬ 
num  roof;  6  =  high-slope,  cedar  rooO 
k  thermal  conductivity  of  plywood  sheathing 

k'  thermal  conductivity  of  shingles 

L  length 

p  angle  (in  horizontal  plane  measured  CW)  between 
ridge  line  and  north 
S'  solar  absorptivity  of  shingles 

W  width 

o  slope  =  low,  high  (for  16.3®,  39.8°) 
r  thickness  of  plywood  sheathing 

r'  thickness  of  shingles  at  butt  end 

Roof  responses 

(rjj;i  =:  E,L,M,R,U,F\  temperature  of  the  /th  section 
/  =  1 ,2,3,4, 5, 6 1  of  the  /th  roof 


Description  of  roofs 

Two  contiguous  gable  roofs  were  built  on  top  of 
a  trailer  being  used  as  an  office  at  CRREL,  Hanover, 
New  Hampshire.  These  roofs  differed  in  slope  but 
were  otherwise  constructed  the  same.  The  sheathing 
of  the  combined  roofs  was  divided  into  seven  identi¬ 
cal  panels  —  three  of  them  on  the  high-slope  end  and 
four  on  the  low-slope  end.  Six  of  these  panels  were 
covered  with  one  of  three  different  roofing  materials  — 
cedar  shingles,  corrugated  aluminum  sheeting,  and 
asphalt  shingles  so  that  six  “test  roofs”  were  created. 
(Hereafter,  “roof’  will  refer  to  one  of  these  six  test 
roofs.)  The  seventh  panel  was  provided  as  a  separa¬ 
tion  to  prevent  the  high-slope  roofs  from  shadowing 
the  low-slope  roofs.  Figure  1  shows  the  test  roofs 
shortly  after  they  were  constructed.  Table  II  gives 
their  characteristics.  The  trailer  was  positioned  as 
shown  in  Figure  2,  which  also  gives  the  location  of 
the  monitoring  stations  for  the  meteorological  data. 
The  floor  of  the  attic  (originally,  the  trailer’s  roof) 
was  insulated  with  Styrofoam  panels  and  access  to 
the  attic  was  provided  by  doors  built  into  the  gables. 


Table  11.  Giaracteiistics  of  roofs. 


a 

T 

k* 

A' 

L 

w 

p 

T 

k 

A 

(°> 

c 

(cm) 

(W/m-K) 

,.,.^lsXW'^) 

S‘ 

(mj 

(m) 

C) 

(cm) 

(W/mK) 

(m^  !$)(.)  or 

16.3 

Asphalt 

0.635 

0.145* 

1.475^ 

0.82  to  0.91^ 

2.438 

l.i29 

43 

1.905* 

0.138* 

1.59* 

AJumiDum 

0.0381 

1 12  to  229** 

480  to  946'* 

0.26* 

2.438 

1.829 

43 

1.905* 

0.138* 

1.59* 

Cedar 

1.27 

0.121^ 

1.826*’^ 

0.9'* 

2.438 

1.829 

43 

1.905* 

0.138* 

1.59* 

39.8 

Asphalt 

0.635 

0.145* 

1.475^ 

0.82  to  0.91* 

2.438 

1.829 

43 

1.905* 

0.138* 

1.59* 

Aluminum 

0.0381 

1 12  to  229^ 

480  to  946“* 

0.26* 

2.438 

1.829 

43 

1.905* 

0.138* 

1.59* 

Cedar 

1.27 

0.121^ 

1.826^'^ 

0.9'* 

2.438 

1.829 

43 

1.905* 

0.138* 

1.59* 

Definitions:  A 
A" 
c 
k 
k* 
L 
P 
S' 
w 

a 

T 

T 


thermal  diffusivity  of  sheathing  (plywood) 

thermal  diffusWity  of  covering 

material  used  to  cover  roof 

thermal  conductivity  of  sheathing  (plywood) 

thermal  conductivity  of  covering 

length  of  roof 

angle  (measured  CW  in  horizontal  plane)  between  ridge  line  of  roof  and  north 
solar  absorptivity  of  covering 
width  of  roof 

slope  (16.3°  =  3.5  in  12;  39.8°  =  10  in  12) 
thickness  of  sheathing  (plywood) 
thickness  of  covering  (butt-end  of  shingles) 


Sources: 

1.  Burch,  D.M.,  B.A.  Peavy  and  F.J.  Powell  (1975)  Comparison  between  measured  computer-predicted  hourly  heating 

and  cooling  energy  requirements  for  an  instrumented  wood-framed  townhouia  subjected  to  laboratory  tests. 
ASHRAE  Trans.,  vol.  81.  part  II,  no.  2  363,  p.  70-88. 

2.  Same  as  (1)  but  assuming  cedar  and  redwood  are  similar. 

3.  Calculated  using  values  in  (1)  for  thermal  conductivity,  specific  heat,  and  density. 

4.  Eckert,  E.R.G.  and  R.M.  Drake,  Jr.  (1950)  Introduction  to  the  transfer  of  heat  and  mass.  New  York,  Toronto, 

London:  McGraw-Hill  Book  Co.,  Inc.,  1st  ed.,  p.  266. 

5.  Strock,  C.  and  R.L.  Koral  (Eds.)  (1965)  Handbook  of  air  conditioning,  heating  and  ventilating.  New  York: 

Industrial  Press,  2nd  ed. 

6.  Brown,  A. 1.  and  S.M.  Marco  (1951) /nrroducrion  fo  hear  rrons/er.  New  York,  Toronto,  London:  McGraw-Hill 

Book  Co.,  Inc.,  2nd  ed.,  p.  68. 

7.  ASHRA '  handbook  of  fundamentals  (1974)  American  Society  of  Heating,  Refrigerating,  and  Air-Conditioning 

Engineers,  Inc.,  New  York,  p.  572.  (Assuming  that  white  oak  is  representative  of  other  woods  and  that  wood 
is  approximately  a  black  body  so  that  emissivity  is  approximately  equal  to  absorptivity.) 


Each  test  roof  was  instrumented  as  follows  (see 
Fig.  3);  A  calibrated  gutter  was  provided  to  measure 
the  volume  of  meltwater  runoff.  Four  thermocouples 
were  installed  in  the  covering  material  on  the  south¬ 
eastern  side  to  measure  the  temperatures  of  the  eaves 
T^,  the  lower  edge  of  the  roof  just  above  the  eaves 
the  middle  of  the  roof  T’m  ,  and  the  ridge  of  the 
roof  Tf^ .  Another  thermocouple  was  installed  under 
the  deck  of  each  roof,  directly  beneath  the  thermo¬ 
couple  measuring  to  measure  the  under-deck 
temperature  T^.  Two  heat  flow  meters  were  in¬ 
stalled,  immediately  adjacent  to  the  thermocouples 
measuring  and  Tm  ,  and  these  measured  the  heat 
flow  through  the  eaves  and  the  heat  flow  through 
the  middle  of  the  roof  respectively. 

Two  additional  thermocouples  were  installed. 
These  were  placed  immediately  above  the  insulated 
attic  floor  so  that  one  was  centered  in  the  attic  space 


under  the  high-slope  test  roofs  and  the  other  was 
centered  under  the  low-slope  test  roofs;  these  thermo¬ 
couples  measured  the  temperature  of  the  attic  floor 
Tp  for  the  two  different  slopes. 

The  cables  for  the  thermocouples  and  heat  flow 
meters  were  strung  from  each  test  roof  to  a  data 
logger  located  in  an  adjacent  building. 

Meteorological  data 

Before  1  November  1972,  meteorological  data 
were  taken  as  follows.  Solar  radiation  was  measured 
with  a  solar  radiometer  that  was  mounted  on  the 
building  housing  the  data  logger.  The  output  of  the 
radiometer  was  coupled  to  a  strip  chart  recorder, 
which  plotted  a  continuous  radiation  curve  as  a  func¬ 
tion  of  time.  These  curves  were  integrated  by  hand  to 
obtain  values  of  (total  radiation  energy  incident 
in  the  vertical  direction)  for  each  day.  Other 


3 


pi 


4 


•  Radiation  Gage 


Table  III.  Data  for  test  days  of  minimuin  windspeed. 


Vertical  Eppley 

Windspeed  V  (kmjh)  Air  temperature  (°C)  Snowfalls  solar  radiation 


Date 

Min  Max  Avg 

Min 

Max 

Avg 

(cm) 

(MJ/m^) 

23  D*c  72 

0 

3.22  1.945 

-  6.67  - 

2.78  - 

4.10 

16.76 

1.46 

24  Dec  72 

0 

6.44  2.28 

-  3.61  - 

1.39  - 

2.78 

16.76 

2.845 

2S  Dec  72 

0 

6.44  3.1  S 

-  1.67  + 

2.22  - 

0.24 

16.76 

0.96 

9  Jan 

73 

0 

4.83  2.01 

-28.89  - 

10.00  - 

20.545 

0 

8.20 

13  Jan 

73 

0 

3.22  1.94 

-21.11  - 

4.17  - 

13.46 

0 

9.83 

14  Jan 

73 

0 

8.05  2.415 

-17.22 

0.00  - 

6.39 

0 

4.895 

IS  Jan 

73 

0 

4.83  2.49 

-  4.72  _+ 

4.17  - 

0.60 

0 

5.48 

Snow  depth  of  /th  roof  (cm) 

Date  h(i  =  1) 

h(i=  2) 

h(i^3) 

h<i“4) 

H(i‘5) 

h(l~6) 

23  Dec  72  1.3 

1.3 

1.3 

1.3 

1.3 

24  Dec  72  1.3 

1.3 

1.3 

1.3 

1.3 

25  Dec  72  1.3 

1.3 

1.3 

1.3 

1.3 

9  Jan  73  21.6 

23.5 

30.5 

14.6 

0 

14.6 

13  Jan  73  21.6 

23.5 

30.5 

13.0 

0 

14.0 

14  Jan  73  21.6 

23.5 

30.5 

13.0 

0 

14.0 

15  Jan  73  20.95 

22.9_ 

29.8 

0 

12.7 

Snow  coverage  of  jth  roof  (%} 

Date  ytj  =  1) 

7(1“  2) 

y(i-3) 

y<i“4) 

y(i  =  5) 

II 

23  Dec  72  100 

100 

100 

100 

100 

100 

24  Dec  72  100 

100 

100 

100 

100 

100 

25  Dec  72  100 

100 

100 

100 

100 

100 

9  Jan  73  87.5 

87.5 

87.5 

81 

0 

81 

13  Jan  73  87.5 

87.5 

87.5 

81 

0 

81 

14  Jan  73  87.5 

87.5 

87.5 

81 

0 

81 

IS  Jan  73  87.5 

87.5 

87.5 

81 

0 

81 

meteorological  data  were  taken  from  the  records  kept 
at  Lebanon  Airport,  West  Lebanon,  N.H.  (13  km  to 
the  south),  and  Dartmouth  College,  Hanover,  N.H. 

(1.6  km  to  the  south). 

After  1  November  1972,  hourly,  on-site  data  for 
all  meteorological  variables  were  provided  by  a 
meteorological  team  attached  to  the  U.S.  Army 
Electronics  Command  and  stationed  at  USACRREL, 
Hanover,  N.H. 

Procedure 

Data  were  collected  from  3  December  1971  to 
20  March  1972  and  from  28  November  1972  to  30 
March  1973.  Thermocouple  and  heat  flow  meter 
data  were  automatically  recorded  on  an  hourly  basis. 
Meteorological  data  were  recorded  on  a  daily  basis 
before  1  November  1972  and,  after  that  date,  on  an 
hourly  basis.  Meltwater  runoff  volumes  were  re¬ 
corded  as  necessary  to  prevent  overflow  of  the  gutter 
Snow  depth  measurements  were  made  daily  on  each 
roof  immediately  following  each  snowfall.  Color 


sUdes  (3S4nm  camera,  Kodachrome  film)  were  taken 
at  regular  intervals  to  record  the  coverage  y,  the 
degree  of  icing  /,  and  the  chronological  history  of 
the  snow  cover  for  each  roof  from  snowfall  to  snow¬ 
fall. 


ANALYSIS 

General 

The  approach  was  to  consider  the  meteorological 
variables  as  interacting  with  the  roof  characteristics 
to  produce  the  roof  responses.  From  this  point  of 
view,  it  was  necessary  to  attempt  to  “separate  vari¬ 
ables.”  The  meteorological  data  were  sorted  to  select 
seven  test  days  for  which  the  windspeed  v  was  the 
lowest  of  all  the  test  days.  This  was  assumed  to  mini¬ 
mize  the  wind  effects  on  the  selected  days,  thus 
allowing  a  study  of  the  effects  of  solar  radiation.  The 
values  of  the  meteorological  variables  for  the  selected 
test  days  are  given  in  Table  III. 
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Table  IV.  Frequency  of  tempcratiui.  profiles:  In  each  profile  expression,  temperature  increases 
from  left  to  right  and  indices  in  parentheses  indicate  equal  temperatures. 


1 

(L.ASP) 

:;o.  of 

test 

Profile 

hours 

E(LMR)UF 

UT 

ER(LM)UF 

30 

ERMLUF 

27 

(ELKR)UF 

23 

E(LR)MUF 

7 

M(LR)UFE 

5 

(LMR)UFE 

3 

(LM)RUFE 

3 

R(LM)UFE 

3 

(ER)(LM)UF 

2 

L(MR)UFE 

2 

MLRUFE 

2 

RMLUFE 

2 

RCLM)EUF 

2 

(EL)(MR)UF 

1 

(ER)LMUF 

1 

EL(MR)UF 

1 

(LR)MUFE 

1 

(LR)MEUF 

1 

LMRUFE 

1 

LMREUF 

1 

(MR)LEUF 

1 

(ME)L(EU)F 

1 

R{LM)UEF 

1 

E(LR)UMF  30 
ELRWIF  25 
(ELU)MRF  20 
E{LU)MRF  i  18 
E{LU)(MR)F  1  lU 
FMULRE  10 
(EL)UfffiF  I  7 
E(LRU)MF  6 
EL(RU)MF  6 
EL(MU)RF  3 
E(LU)RMF  2 
(ELU)(MR)F  2 
ELU(MR)F  2 
(MF)ULRE  2 
(MF)(LU)(ER)|  2 
(EF)(LU)MF  1 
(EL)U(MR)F  1 


ERLMUF 

(EL)(MR)UF 

E(I>IR)UF 

E(LR)MUF 

E(LM)(RU)F 

ER(LM(UF 

ELRMUF 

EL(MR)UF 

E(LM)URF 

(ELR)MUF 

LMRUFE 

(ELM)URF 

LMREUF 

LMHUEF 

E(LM)RUF 

(LR)MEUF 

LEMRUF 

(ER)(LM)UF 

ERMLUF 

ELMRUF 

LMR(EU)F 

L(MR)(EU)F 

L(MR)UFE 

LMERUF 

L(EM)RUF 

R(LM)UEF 

RELMUF 

R(EL)MirF 


(ELMU)(RF) 

(EL)(MU)RF 

ERLUMF 

EL(MRU)F 

(LR)UEMF 

(LR)(EMU)F 

MU(ELF)R 

M(EL)UFB 

MF(ELU)R 

(MUF)(LR)E 

MFULER 

(MF)(LU)RE 

R(EL)UMF 

FMUL(ER) 

F(MU)RLE 

FM(LU)(ER) 

FM(LU)RE 


EL(MR)UF  68 
ELRMFU  30 
EL(MR)FU  27 
ELRFMU  8 
ELR(MF)U  6 
RMLEFU  5 
RMLFEU  k 
E(LMR)UF  3 
EL(MR)(UF)  3 
RMLFUE  3 
(ER)LMFU  2 
ELRUMF  1 
(ER)(LM)FU  1 
ELMR(UF)  1 
(ELR)MFU  1 
(EL)(MR)FU  1 
(EL)(MR)(UF)  1 
(MR){EL)FU  1 
R(LM)EFU  1 
RfELM)FU  1 


ELM(RU)F 

ERLMUF 

ELMRUF 

ER(LM)UF 

MLUFRE 

EL(MR)UF 

LMRUEF 

ERMULF 

(LM)URFE 

(lmHer)uf 

(LMR)EUF 

ML(ERU)F 

R(MU)LFE 

(EL)MRUF 

ELMURF 

(ER)MU(LF) 

ER(MU)(LF) 

(ER)(LM)UF 

E(LR)MUF 

(ELM)RUF 


(H,ASP) 

IJo.  of 

test 

Profile 

hours 

E(LM)UF 

ELMUF 

(ELM)UF 

(MUFLE 

MLUFE 

M(UF)LE 

(EL)MUF 

EL(MUF) 

EL(MU)F 

MFULE 

MU(LF)E 

MULFE 

ELU(MF) 

EM(UF)L 

EMU(LF) 

LM(EU)F 

LMUFE 

(LM)UEF 

MLEUF 

M(LU)FE 

MLU(EF) 


No.  of 
test 

Profile  hours 


ERMLUF 

L(MR)UEF 

(LM)RUFE 

LM(RU)EF 

(LM)ERUF 

L(EM)RUF 

(LM)(RU)EF 

(LM)UREF 

M(LRU)FE 

MLU(RF)E 

MULRFE 

(MR)ULFE 

MLU(ER)F 

R(MU)E(LF) 

RE(MU)(LF) 

RMULFE 

RM(LU)EF 

REMULF 

R(ELM)UF 


Table  V.  Maximum  temperature  profiles.  In  each  profile  expression,  temperature  increases 
from  left  to  right. 


Roof,  j 

Maximum 

Profile 

No.  of 
test 
hours 

Observed  Members  of  Compatible  Set 

1 

ERMLUF 

129 

E(LMR)UF,  ER(LM)UF,  ERMLUF,  (ELMR)UF, 
(ER)(LM)UF 

2 

ELRMUF* 

T9 

(EL)(MR)UF,  E(LMR)UF,  E  (LR)MUF, 

ELRMUF,  EL(MR)UF,  (ELR)MUF 

ERLMUF 

T1 

ERLMUF,  E(LMR)UF,  E(LR)MUF,  ER(LM)UF, 
(ELR)MUF,  (ER)(IM)UF 

3 

ELMRUF* 

76 

EL(MR)UF,  E(LMR)UF,  EL(MR)(UF), 

ELMR(UF),  (EL)(MR)(UF) 

ELRMUF 

75 

EL(MR)UF,  E(LMR)UF,  EL(MR)(UF), 
(EL)(MR)(UF) 

ELRMFU 

69 

ELRMFU,  EL(MR)FU,  ELR(MF)U,  EL(MR)(UF), 
(ELR)MFU,  (EL)(MR)FU,  {EL)(MR)(UF) 

k 

ELMUF 

137 

E(LM)UF,  ELMUF,  (ELM)UF,  (EL)MUF, 
el(muf),  el(mu)f 

5 

ELUMRF* 

70 

(ELU)MRF,  E(LU)MRF,  E(LU)(MR)F, 

(el)umrf,  el(mu)rf,  {elu)(mr)f, 
elu(mr)f,  (el)u(mr)f,  (elmu)(rf), 
(el)(mu)rf,  el(mru)f 

ELRUMF 

68 

E(LR)UMF,  ELRUMF,  E(LRU)MF,  EL(RU)MF, 
EL(MRU)F 

6 

ELMRUF 

83 

ELM(RU)F,  ELMRUF,  EL(MR)UF,  (EL)MRUF, 
(ELM)RUF 

•Significant  observations  of  competing  profiles  make  this  identification 
as  "maximum"  unsure.  In  these  oases,  the  ccanpeting  profiles  are  given 
below  the  tentative  maximum  profile. 


Temperature  profiles 

For  a  given  roof  /  and  a  given  hour  of  the  day  t,  the 
temperatures  r^j,  T^jj,  and  Tf,  (where  the 

indices  refer  to  the  roof  section  and  roof  as  defined 
in  Table  I)  were  compared  to  provide  a  “profile”  of 
the  roof  for  that  hour.  The  profile  was  symbolized 
by  writing  the  indices E,  L,M,R,  U,  and  F  in  the 
order  of  increasing  temperature  values  with  equal 
temperatures  represented  by  parentheses  around  the 
relevant  indices.  For  example,  “FL(URfME"  meant 
that,  on  a  given  day  at  a  given  hour  for  roof/,  Tj.,  < 
^Lj  <  ^UJ  “  ^Rj  <  ^MJ  <  ^Ejv  profiles  were 
and  as  convenient  shorthand  to  display  the  manner 
in  vdiich  the  different  parts  of  a  roof  changed  in  temp¬ 
erature  relative  to  each  other  as  the  solar  radiation 
changed.  Profiles  were  written  for  each  roof  and 
each  of  the  168  hours  of  the  seven  test  days.  For 
each  roof,  the  number  of  hours  for  which  each  profile 
was  observed  was  counted.  Table  IV  lists  the  profiles 
observed  and  their  frequency  for  each  roof. 


Once  the  distribution  of  profiles  had  been  obtained, 
this  distribution  was  analyzed  to  determine  the  profile 
containing  no  equalities  that  most  nearly  accounted  for 
all  the  profiles  observed  for  a  given  roof.  For  example. 
Table  IV  shows  that  the  profiles  (ELYMR)UF, 
E(Jjm)VF,  E{LR)MUF,  E(Lhr0.U)F,  ELRMUF, 
EL{MR)UF,  and  (ELR)MUFut  all  compatible  in  that, 
for  the  ( qualities,  the  order  of  indices  is  immaterial. 
This  entire  set  of  profiles  can  be  represented  by 
ELRMLF,  the  member  with  no  equalities.  All  of  the 
possible  compatible  sets  of  this  sort  were  assembled 
for  each  roof  and  the  total  numbers  of  observations 
fsr  each  counted.  (To  distinguish  these  sets  of  pro¬ 
files  from  the  profile  that  is  only  a  member  of  the  set, 
the  profile  expression  was  underlined  when  it  re¬ 
ferred  to  the  compatible  set.)  The  compatible  set 
with  the  largest  number  of  hours  of  observation  was 
termed  the  “maximum  profile”  for  the  roof.  Table  V 
lists  the  maximiun  profiles  for  the  six  test  roofs. 
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Tabte  VI.  Temperature  profile  variation  with  solar  radiation. 

The  “most  likely’*  profile  for  each  hour  of  a  day  it  compared  with  the  values  of  the  hourly  tolar  radiation  for 
14  January  1973.  (The  total  daily  solar  radiation  was  a  median  value  for  14  January  1973  as  compared  with 
the  other  six  test  days.)  In  each  profile  expression,  temperatures  increased  from  left  to  right. 


MOST 

LIKELY  PROFILES  FOR 

jth  ROOF 

1 

2 

3 

1* 

5 

6 

n 

Hour 

(L,ASP) 

{L,ALM) 

(L,CDR) 

(H,ASP) 

(H,ALM) 

(H,CDR) 

MJ/m^ 

0100 

ERMLUF 

ERLMUF 

ELRMFU 

ELMUF 

ELRUMF 

ERLMUF* 

0 

0200 

It 

ELRMUF 

«  * 

" 

It 

tl 

tl 

0300 

" 

ERLMUF 

»•  » 

II 

II 

tl 

tl 

OhOO 

tl 

II 

II 

II 

II 

tl 

0500 

II 

It 

»'  « 

It 

tt 

It 

’’ 

o6oo 

II 

II 

"  II 

It 

It 

II 

tl 

0700 

*' 

M 

»  H 

" 

tl 

It 

II 

o8oo 

tl 

ELRMUF* 

*•  # 

It 

II 

tl 

0.25 

0900 

II 

»*  « 

It 

II 

ELUMRF 

ELMRUF 

0.29 

1000 

II 

II 

tt 

II 

tl 

ft 

0.51* 

1100 

»»  « 

ELMURF 

ELMRUF 

♦'  » 

tl 

It 

0.83 

1200 

”  * 

It 

It 

MUFLE 

FMULRE 

II 

0.96 

1300 

'»  # 

tl 

'*  H 

tl 

It 

0.63 

litOO 

"  * 

II 

"  n 

11 

It 

0.5I* 

1500 

II 

It 

It 

ELMUF 

It 

II 

0.50 

i6oo 

II 

" 

It 

It 

ELUMRF 

" 

0.25 

1700 

It 

ELRMUF 

If 

If 

If 

M 

0.08 

1800 

II 

II 

tl 

tl 

»•  * 

It 

0 

1900 

II 

ERLMUF* 

tl 

tl 

ELRUMF 

II 

II 

2000 

tl 

'»  * 

tt 

II 

tl 

II 

2100 

" 

It 

ELRMFU 

tt 

If 

tl 

,, 

2200 

M 

ELRMUF 

tt 

It 

tl 

II 

tl 

2300 

If 

tl 

tl 

It 

11 

tt 

It 

21*00 

tl 

ERLMUF* 

tl 

II 

It 

It 

II 

•  Significant  observations  of  competing  profiles  make  this 
identification  as*)nost  likely^^unsure. 


The  profiles  were  then  considered  by  time  of  day 
for  each  roof.  For  each  hour  of  the  day  and  a  given 
roof,  a  profile  for  each  of  the  seven  test  days  was 
written.  The  compatible  set  to  which  most  of  these 
seven  profiles  belonged  was  then  termed  the  “most 
likely  profile”  for  that  roof  at  that  time  of  day.  For 
example,  if  the  low-slope  aluminum  roof  exhibited  the 
following  profiles  at  0800  hours  on  the  seven  test 
days:  {ELR)MUF,  ELRMUF,  ERLMUF,  ERLMUF, 
ELRMUF,  ELRMUF,  and  {EL)iMR)UF,  the  most 
likely,  profile  for  this  roof  at  0800  hours  was  ELRMUF. 
The  results  of  this  analysis  are  given  in  Table  VI,  which 
displays  the  most  likely  profiles  as  a  function  of  roof, 
time  of  day,  and  solar  radiation. 


Temperature  rankings 

For  a  given  section  i  of  a  roof  and  a  given  hour  of 
the  day,  the  temperatures  Tii,  712,  ri3.  ri4,  7)$,  and 
Tie  (where  the  numerical  indices  now  refer  to  the  roof 
number  /  as  defined  in  Table  I)  were  compared  to  pro¬ 
vide  a  “ranking”  of  the  roofs  for  that  section  at  that 
hour  of  that  day.  These  rankings  were  expressed  in  a 
manner  analogous  to  that  discussed  above  for  profiles 
and  analyzed  in  the  same  manner  to  find  the  maximum 
ranking  for  a  roof  part  and  the  most  likely  ranking  for  a 
roof  part  at  a  given  time  of  day.  For  example,  the  rank¬ 
ing  126543  meant  that  Tjj  <  Tjj  <  <  Tjj  <  Tj^  < 

T’.s  and  the  ranking  1 26543  represented  a  set  of  rank¬ 
ings  compatible  with  126543. 


Table  Vll.  Frequency  of  temperature  ranking?  for  tth  roof  section. 

The  lists  of  rankings  have  been  abbreviated  to  give  only  the  first  16  rankings  in  order  of  decreasing  frequency. 
This  is  sufTicient  to  account  for  at  least  S0%  of  the  data  hours. 


E  =  EAVES 

L  =  LOW 

H  =  MIDDLE 

R  =  RIDGE 

Ranking 

No. 

Ranking 

No. 

Ranking 

No. 

Rank! ng 

No. 

(65‘<321) 

16 

(651*321) 

23 

(1*21)3(65) 

56 

13265 

23 

3(651*21) 

9 

51*63(21) 

13 

(1*21)365 

7 

65123 

21 

(61*321)5 

8 

2(651*31) 

13 

(61*)1(52)3 

5 

(21)356 

15 

326511* 

T 

(61*321)5 

10 

(61*)5123 

5 

1(32)65 

9 

(632)(5l*l) 

5 

(21)631*5 

8 

61*(51)23 

5 

56123 

8 

362511* 

5 

3(651*21) 

5 

5(61*)123 

5 

(65)123 

6 

(6l4  32)(5l) 

1* 

21631*5 

I4 

(1.21)(53)6 

5 

6(51)23 

6 

(63)(5l*2l) 

1* 

(32)(61*51) 

1* 

6(1*1)523 

U 

213(65) 

6 

(632)(1.1)5 

1* 

(51*321)6 

1* 

6(1*1)2(53) 

3 

(21)3(65) 

6 

1*15263 

1* 

1*(653)(21) 

1* 

(61)21*35 

3 

132(65) 

5 

1*51236 

3 

32(651*1) 

3 

6(1*1)(52)3 

3 

123(65) 

5 

3(61)(51*2) 

3 

(51*3)621 

3 

6(U21)35 

3 

12365 

1* 

31(651*2) 

3 

(51*3)6(21) 

3 

6U1253 

2 

16235 

1. 

36(51*21) 

3 

51*6312 

3 

(61*)1253 

2 

(65)(21)3 

3 

3(61*21)5 

3 

51*(21)36 

3 

6(51*1)23 

2 

65(21)3 

3 

3(62)511* 

3 

(62)(1*31)5 

3 

6(1*1)235 

2 

(6i)(52)3 

3 

others 

81* 

others 

62 

others 

56 

others 

1*1 

tJ  =  UNDER 

DECK 

F  =  FLOOR 

Ranking 

No. 

Ranking 

No. 

561*123 

1*3 

(65l*)(32l) 

95 

52l*(6l)3 

11 

(651*321) 

57 

561*(21)3 

11 

(321)(651*) 

16 

5(l*2)(6l)3 

8 

5l*(21)63 

8 

51*(621)3 

7 

5l*(21)(63) 

7 

(65)1*123 

6 

52(1*1)63 

5 

5(1*2)163 

5 

61(1*2)35 

1* 

61*1253 

1* 

(5l*2)(6l)3 

1* 

61*1235 

3 

5(1*2)1(63) 

3 

51*21(63) 

3 

others 

36 

It  should  be  noted  at  this  point  that,  for  the  floor  of 
the  roof,  only  two  thermocouples  were  installed  so 
that  the  high-slope  and  the  low-slope  roofs  could  be 
distinguished  but  not  the  individual  high-slope  roofs  or 
the  individual  low-slope  roofs.  That  is,  for  the  floor 
(i  =  F),  there  were  only  three  possible  rankings:  (321) 
(654),  (654X321),  and  (654321). 

Table  VII  gives  an  abbreviated  listing  of  the  ob¬ 
served  rankings  for  each  roof  section.  Table  VIII  lists 
the  maximum  rankings  for  each  roof  section  and  the 
most  likely  rankings  are  given  in  Table  IX  as  a  function 
of  the  hour  of  the  day  and  the  solar  radiation. 

A  comparison  of  Tables  IV,  V  and  VI  with  Tables  VII, 
VIII,  and  IX  immediately  indicates  that  the  interpretation 
of  the  variation  for  a  given  roof  is  more  straightforward 

9 


than  that  for  the  variation  for  a  given  roof  section.  This 
is  because,  in  the  latter  case,  there  are  competing  effects 
of  roof  covering  composition  c  and  roof  slope  a.  The 
effects  of  these  two  roof  characteristics  must  be  separated. 

To  determine  the  relative  effects  of  the  roof  character¬ 
istics  a  and  c\  Tables  X  and  XI  were  constructed.  Table 
X  gives  the  percentage  of  the  total  number  of  test  hours 
for  which  each  of  the  possible  ranking  equalities  was 
observed  for  a  given  roof  section.  This  displays  which 
roofs  tend  to  behave  alike  for  a  given  roof  section  and 
should  reveal  effects  due  to  composition.  Table  XI 
was  constructed  by  first  eliminating  all  those  rankings 
for  which  the  presence  of  equalities  made  it  unclear 
whether  the  high-slope  roofs  as  a  group  were  lower  or 
higher  in  temperature  for  the  given  roof  section  than 


i.viw.. 


Table  VIII.  Maximum  temperature  rankings. 


ith 

Section 

MaximuB 

Ranking 

No. 

Observed  Members  of  Compatible  Set 

E  =  Eaves 

3621*15 

61 

(651*321).  3(651*21),  (61*321)5,  (632)(51*1) 
(61*32)(51),  (63)(5l*2l),  (632)(i*l)5, 
36(1*21)5,  3(62)(5l*l),  3(61*2)(51) 

L  =  Low 

21631*5 

62 

(651*321),  2(651*31),  (61*321)5,  (21)631*5, 
21631*5,  (621)31*5,  2163(51*),  2(61*31)5 

M  =  Middle 

61* 

(1*21)3(65),  (1*21)365,  (1*21)(63)5 

R  =  Ridge* 

65123 

1*0 

65123,  (65)123,  6(51)23,  (65)(21)3, 
65(21)3,  (651)23 

13265 

39 

13265,  1(32)65,  132(65),  1(632)5 

U  =  Under 
deck 

561*123 

63 

561*123,  561*(21)3,  (65)1*123,  56(1*1)23 

95 

(65l*)(321) 

*More  than  one  set  included  due  to  lack  of  distinction  in  number  of  observed  hours  of  data. 


the  low-slope  roofs.  The  remaining  rankings  were  di¬ 
vided  into  three  groups:  those  with  all  high-slope 
roofs  having  higher  temperatures  than  all  the  low-slope 
roofs,  those  for  the  reverse  case,  and  those  for  which 
neither  case  was  true.  This  displays  effects  due  to  the 
slope  of  the  roof. 


time  periods  f ,  to  f  j ,  0000  to  0800  hours,  0800  to 
1700  hours,  and  1700  to  2400  hours,  of  two  repre¬ 
sentative  test  days.  The  results  are  given  in  Table  XII 
along  with  the  relevant  meteorological  data  and  the 
corresponding  average  temperatures  of  the  eaves  and 
middle  roof  sections. 


Heat  flow 

The  net  heat  accumulated  per  unit  area  in  the  time 
period  f  j  to  fj  by  the  eave  (or  middle)  sections  was 
calculated  from  the  heat  flow  meter  data  for  the  seven 
test  days  by  use  of  the  approximation: 


where  is  the  value  of  that  was  measured  at 
the  beginning  of  hour  t.  (The  equivalent  expression 
for  the  middle  section  is  achieved  by  substituting 
and  for  and ^^(r).]  This  approximation 
means  to  assume  that  the  heat  flow  (or  re¬ 
mained  constant  during  any  given  test  hour,  ^uation 
1  is  written  so  that  net  heat  flow  into  the  eaves  (or 
middle)  section  of  a  roof  will  be  indicated  by  a  posi¬ 
tive  value  {otJUy^dt).  The  net  heat  flows 

Sn^dt  vnAfiiy^dt  were  calculated  for  three  speciflc 


Snow  depths  and  meltwater  volumes 

The  snow  depths  Aj,  Aj,  Aj,  A^,  Aj  and  A^  (where 
the  indices  refer  to  the  particular  test  rooO  that  were 
observed  on  the  roofs  are  given  in  Table  XIII  and  the 
corresponding  percentages  of  cover  7^,  72,  73,  7^, 

7; ,  and  7*  are  given  in  Table  XIV. 

The  fraction  of  the  original  snow  cover  that  was 
collected  as  meltwater  volumes  K,,  Ij,  Fj,  F4,  Kj 
and  where  the  indices  again  refer  to  the  test  roofs 
was  crJculated  by  assuming  that  the  original  snow  cover 
had  a  density  of  0.2  g/cm*  and  had  a  unifonn  thick¬ 
ness.*  “Original  snowfall”  is  the  snowfall  occurring 
in  such  close  proximity  in  time  as  to  be  undissipated 
at  the  time  of  the  first  measurements  of  snow  depths 
and  coverages  following  its  cessation.  The  meltwater 
volumes  were  measured  wiienever  the  practical  re¬ 
quirements  of  using  the  calibrated  gutters  made  it 
necessary.  With  these  assumptions,  the  meltwater 
fraction  of  the  original  snow  cover  could  be  calculated 
for  each  roof  from 


Table  IX.  Temperature  ranking  variation  with  solar  radiation. 

The  “most  likely”  ranking  for  each  hour  of  a  day  is  compared  with  the  values  of  the  hourly  solar  radiation  for 
14  January  1973.  (The  total  daily  solar  radiation  was  a  median  value  for  14  January  1973  as  compared  to 
the  other  six  test  days.) 


Hour 

J 

K- Haven 

MOST  n 

L=Lo« 

KELY  KANKIjffl 

M=Middle 

rrcRTtiT 

R=Rldge 

ROW  Sect] 

U=Under 

roN 

F=Floor 

c 

vu%^ 

0100 

1)51632 

51)6321 

(1)21)3(65) 

65123 

561)123 

(65l))32l) 

0 

0200 

tl 

II 

II 

" 

II 

II 

II 

0300 

1)16352 

II 

II 

•1 

II 

II 

II 

OliOO 

1)15632 

II 

II 

II 

11 

II 

II 

0500 

M 

II 

II 

It 

II 

II 

0600 

161)352 

II 

II 

II 

0700 

II 

II 

II 

It 

II 

0800 

1631)25 

II 

II 

■■ 

II 

0.25 

0900 

3621)15 

tl 

It 

HI 

IB 

II 

0.29 

1000 

36215!) 

II 

II 

Bi 

II 

0.51) 

1100 

tl 

62131)5 

II 

It 

0.83 

1200 

II 

II 

II 

13265 

61)1235 

II 

0.96 

1300 

II 

2l63l)5 

II 

" 

6ll)235 

It 

0.63 

1^*00 

" 

■1 

It 

" 

II 

"  « 

0.51) 

1500 

32651I) 

mm 

II 

II 

II 

0.50 

1600 

3261)15 

tl 

»'  * 

521)163 

0.25 

1700 

321)165 

II 

II 

12365 

»*  # 

0.08 

1800 

321)615 

II 

II 

561)123 

0 

1900 

3261)15 

■1 

II 

" 

" 

It 

2000 

1)12536 

rt 

" 

II 

II 

2100 

3261)15 

mm 

•1 

65123 

II 

M 

2200 

1)15263 

B 

II 

II 

"  # 

2300 

11)5263 

II 

" 

" 

" 

2U00 

1)15623 

II 

II 

II 

II 

_ ^ _ 

It 

_ 

•  Significant  observations  of  competing  rankings  make  this  identification 
as  "most  likely"  unsure. 


Meltwater 
Snow  cover 


roof  no. 

(2) 


and  the  results  are  listed  along  with  the  relevant  melt¬ 
water  volumes  in  Table  XV.  In  eq  2,  L  and  W  are  the 
common  length  and  width  of  the  test  roofs  (refer  to 
Table  10;  the  other  variables  are  meltwater  volume  V, 
coverage  y,  and  snow-depth  h,  which  arc  functions  of 
the  particular  roof  and  meteorological  variables. 

In  Table  XV,  all  data  are  keyed  to  the  “original 
snowfall!”  Upon  cessation  of  a  snowfall,  measure¬ 
ments  of  h  and  y  for  each  roof  were  made  at  intervals 


to  provide  a  record  of  the  changing  conditions  on  the 
roofs  following  snowfall.  The  values  of  h  and  y  that 
were  used  in  eq  2  to  compute  meltwater  fractions 
were  those  measured  immediately  after  cessation  of 
the  snowfall  before  significant  dissipation  had  occurred. 
Depending  upon  meteorological  conditions,  the  cali¬ 
brated  gutters  required  emptying  at  differing  fre¬ 
quencies  for  different  snowfalls.  The  meltwater  frac¬ 
tions  and  volumes  recorded  in  Table  XV  represent  the 
total  meltwater  fraction  for  the  particular  snowfall: 

The  relevant  values  of  V  for  each  roof  were  determined 
by  adding  all  of  the  meltwater  volumes  collected  for  a 
given  roof  subsequent  to  the  given  snowfall.  The  dates 
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Table  X.  Effects  due  to  composition  of  roof  covering. 

For  each  tectlon,  the  percenUfe  of  tett  hours  for  which  the  given  roofs  were  equal  in  temperature  is  displayed. 
The  equal  temperature  roofs  are  indicated  by  the  equality  portions  of  rankings.  All  possibilities  for  equalities 
are  listed. 


Partial 

Ranking 

PER 

E= Eaves 

CENT  OF  TE 
L=Low 

ST  HOURS  OB 
M=Middle 

SERVED  FOR  i 
R=Ridge 

th  ROOF  SECT 
U=Under 

ION 

F=Floor 

(651*321) 

9.5 

llt.O 

0 

0 

0 

31*.  0 

(651*32) 

tl 

It 

II 

It 

tl 

It 

(651*31) 

10.0 

21.0 

It 

It 

It 

It 

(651*21) 

15.0 

18.0 

II 

tt 

tl 

fl 

(65321) 

ll*.0 

20.0 

tl 

II 

tl 

tl 

(51*321) 

11.0 

16.0 

tt 

tt 

It 

tt 

(651*3) 

10.0 

23.0 

II 

tl 

It 

II 

(651*2) 

IT.O 

19.0 

It 

II 

II 

tl 

(651*1) 

15.0 

27.0 

0.5 

M 

II 

tl 

(6532) 

11.0 

15.0 

0 

It 

tl 

tl 

(6531) 

10.0 

22.0 

tl 

II 

tl 

tl 

(61*32) 

18.0 

20.0 

It 

It 

tl 

tl 

(61*31) 

15.0 

28.0 

It 

It 

It 

tl 

(51*32) 

11.0 

17.0 

It 

II 

*t 

II 

(51*31) 

M 

26.0 

It 

It 

tl 

tt 

(1*321) 

16.0 

23.0 

tl 

tl 

tt 

II 

(5321) 

11.0 

17.0 

It 

" 

rt 

tt 

(51*21) 

20.0 

21.0 

It 

0.5 

tl 

(6321) 

15.0 

20.0 

tl 

It 

0 

II 

(6U21) 

21.0 

2U.0 

It 

II 

II 

(6521) 

15.0 

19.0 

II 

tl 

II 

It 

(651*) 

17.0 

31.5 

1.0 

tt 

II 

100.0 

(61*3) 

19.0 

30.0 

0 

II 

II 

31*. 0 

(632) 

26.0 

23.0 

It 

1.0 

It 

tl 

(621) 

23.0 

26.0 

It 

It 

5.0 

tt 

(653) 

12.0 

2I1.O 

" 

0 

0 

tt 

(652) 

18.0 

20.0 

II 

0.5 

II 

II 

(651) 

15.0 

31.0 

1.0 

II 

It 

It 

(61*2) 

28.5 

25.0 

tl 

0 

It 

tt 

(61*1) 

23.0 

37.0 

2.0 

tl 

0.5 

It 

(631) 

15.0 

28.5 

0 

It 

ft 

It 

(51*3) 

11.0 

31.5 

It 

tl 

0 

It 

(51*2) 

21*.  0 

2.1 

tt 

It 

5.0 

tl 

(51*1) 

25.0 

35.0 

2.0 

tl 

0 

tl 

(532) 

15.0 

18.0 

0 

tl 

tl 

It 

(531) 

11.0 

27.0 

tl 

tl 

It 

tl 

(521) 

20.0 

21.0 

It 

0.5 

1.0 

tt 

(1*32) 

21.0 

23.0 

II 

0 

0 

It 

(1*31) 

22.0 

37.0 

tl 

11 

tl 

tl 

(1*21) 

29.0 

29.0 

1*5.0 

tl 

1.0 

II 

(321) 

17.0 

25.0 

0.5 

It 

0 

100.0 

(65) 

22.0 

36.0 

37.0 

20.0 

3.5 

100.0 

(61*) 

33.0 

1*2.0 

15.0 

0 

0.5 

tl 

(63) 

tl 

1*0.0 

3.0 

2.0 

11.0 

3l*.0 

(62) 

1*0.0 

31*. 5 

1*.0 

5.0 

5.0 

tl 

(61) 

28.5 

39.0 

8.0 

It 

21.0 

tt 

(51*) 

32.0 

1*6.0 

3.0 

0 

7.0 

100.0 

(53) 

20.0 

37.0 

8.0 

It.o 

0.5 

31*. 0 

(52) 

33.0 

23.0 

6.0 

tl 

7.0 

3l*.0 

(51) 

31.0 

36.0 

8.0 

6.5 

2.0 

tt 

(U3) 

21.0 

1*1*. 0 

2.0 

0 

0 

II 

(U2) 

38.0 

30.0 

1*9.0 

tl 

19.0 

It 

(1*1) 

1*2.0 

1*5.0 

57.0 

It 

10.0 

If 

(32) 

35.0 

25.0 

2.0 

6.5 

0 

100.0 

(31) 

17.0 

1*0.0 

0.5 

0 

0.5 

tt 

(21) 

32.0 

55.0 

1*6.0 

22.0 

2l*.0 

ft 
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Table  XI.  Effects  due  to  slope  of  roof. 

For  each  section,  the  percentage  of  test  hours  for  which  the  high>slope  roofs  as  a  group  were  definitively 
lower  or  higher  in  temperature  than  the  low-slope  roofs  as  a  group  is  tabulated.  All  data  for  which  the  pre* 
sence  of  equal  temperatures  made  it  unclear  whether  the  high-slope  or  low-slope  roofs  were  higher  in  temp¬ 
erature  were  excluded  from  consideration. 


Roof 

PERCENT  OF  TEST  HOURS  OBSERVED 

section 

High-slope  higher 

Low-slope  higher 

Equal 

Eaves 

0 

0 

9.5 

Low 

28.0 

0 

lU.9 

Middle 

11.0 

0 

0 

Ridge 

27.0 

51.0 

0 

Under 

1*2.0 

0 

0 

Floor 

56.5 

9.5 

31*. 0 

given  are  those  spanning  the  time  from  cessation  of  the 
snowfall  to  dissipation  of  the  snow/ice  from  the  roofs 
or  the  beginning  of  a  new  snowfall,  whichever  applied. 
The  average  total  solar  radiation  for  this  period  is  also 
given. 

Degree  of  icing 

A  semi-quantitative  estimate  of  the  degree  of  icing 
/  was  made  by  using  the  definition  of  /  in  Table  1  to 
analyze  the  chronological  sequence  of  photographic 
data.  The  results  of  this  analysis  are  in  Table  XVI. 

The  relevant  photographs  are  illustrated  in  Figure  4. 


RESULTS  AND  CONCLUSIONS 
General 

In  all  discussion  below,  the  test  roofs  will  be  identi¬ 
fied  both  descriptively  and  numerically  as  in  the  ex¬ 
pression  “no.  3  (low-slope,  cedar)  roof.”  For  more  de¬ 
tailed  information,  refer  to  Tables  I  and  II. 

Individual  results  will  be  recounted  followed  by  a 
discussion  relating  them  to  the  icing  problem. 

Temperature  profiles  of  the  roofs 

Analysis  of  the  results  compiled  in  Table  IV  allows 
the  roofs  to  be  ordered  according  to  the  stability  of 
their  temperature  profiles,  i.e.,  their  ability  to  main¬ 
tain  the  same  temperature  profile.  In  order  of  decreas¬ 
ing  stability: 

no.  3  (low-slope,  cedar)  roof 
no.  4  (high-slope,  asphalt)  roof 
no.  6  (high-slope,  cedar)  roof 
no.  1  (low-slope,  asphalt)  roof 


no.  S  (high-slope,  aluminum)  roof 
no.  2  (low-slope,  aluminum)  roof 
With  the  exception  of  the  cedar  roofs,  this  order 
implies  that  higher  slope  leads  to  a  more  stable  temp¬ 
erature  profile.  The  reason  that  the  cedar  roofs  fail  to 
display  the  same  slope  dependence  is  that  the  no.  3 
(low-slope,  cedar)  roof  was  the  only  roof  that  retained 
a  deep  enough  snow  cover  (approaching  30  cm)  to 
screen  the  effects  of  solar  radiation.’  A  study  of  Table 
II  will  reveal  that  the  cedar  shingles  have  the  thickest 
butt  end,  r',  and  therefore,  for  a  given  slope,  the  great¬ 
est  ability  to  retain  snow  cover  from  a  purely  geo¬ 
metrical  view. 

The  stability  of  the  temperature  profiles  as  defined 
above  does  not  take  into  account  how  great  are  the 
changes  in  temperature  profile;  that  is,  how  significant 
are  the  changes  in  profile  irrespective  of  how  often  the 
changes  take  place?  This  question  can  be  answered  by 
using  the  concept  of  compatible  sets  of  profiles.  Those 
profiles  belonging  to  a  compatible  set  are  not  signifi¬ 
cantly  different.  The  maximum  profiles  are  compatible 
sets  most  often  observed.  Table  V  records  the  maxi¬ 
mum  profiles  and  may  be  used  to  obtain  an  ordering 
of  the  roofs  according  to  their  ability  to  maintain  com¬ 
patible  temperature  profiles: 

no.  4  (high-slope,  asphalt)  roof 
no.  1  (low-slope,  asphalt)  roof 
no.  6  (high-slope,  cedar)  roof 
no.  2  (low-slope,  aluminum)  roof 
no.  3  Qow-slope,  cedar)  roof 
no.  5  (high-slope,  aluminum)  roof 
Again,  the  results  imply  that  the  greater  the  slope, 
the  greater  the  ability  to  maintain  a  temperature  pro¬ 
file  within  a  narrow  range  of  change.  This  time,  it  is 
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4a.  0830.  1  Dec  71. 

Figure  4.  Photographs  displaying  relative  icing  of  the  roofs. 
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4g.  1 245.  5  Jan  72. 


Table  XII.  Heat  accumulated  by  eaves  and  middle  roof  sections. 

Data  are  for  three  daily  time  periods  and  two  test  days. 


Date 

Averai?:e 

temperature 

{°C) 

Net  heat  flow 
(MJ/r^) 

Solar 

radiation 

Max . 
wind- 
speed 

Avp. 

ai  r 
temp. 

Time 

Roof 

Eaves 

Middle 

Eaves 

Middle 

(MJ/n^) 

(kiti/h ) 

(‘"c) 

23 

1 

-1  J,1 

+0.00 

-0.028 

-o.oi*i* 

0 

3.2 

-5.7 

Dec 

T2 

2 

-0.555 

II 

+0.007 

+0.133 

tl 

II 

" 

0000 

3 

-1.91) 

+0.555 

+0.001 

+0.162 

ri 

II 

II 

to 

0800 

-0.555 

+0.00 

+0.032 

+0.11*1 

II 

II 

II 

5 

.. 

+0.555 

-O.OU5 

+0.162 

II 

II 

6 

-0.83 

+1.11 

+0.019 

+0.11*3 

11 

II 

II 

0800 

1 

-0.555 

+0.00 

-0.019 

+0.02I* 

1.1*6 

-3.55 

to 

1700 

2 

rt 

II 

+0.126 

+0.21*1 

II 

If 

•1 

3 

-1.67 

+0.555 

-0.098 

+0.2725 

II 

II 

1* 

-0.555 

+0.00 

+0.112 

+0.190 

II 

II 

l« 

5 

+0.00 

+1.11 

-0.018 

+0.187 

II 

II 

II 

6 

-0.555 

II 

+0.101 

+0.2225 

II 

II 

II 

1700 

1 

It 

-0.555 

-0.021t 

-0.052 

0 

-3.3 

to 

2^*00 

2 

M 

+0.00 

+0.056 

+0.11*8 

II 

3 

-1.11 

+0.555 

-0.265 

+0.177 

II 

1. 

-0.555 

+0.00 

+0.0it8 

+0.151 

5 

It 

+1.11 

-0 . 020 

+0.178 

II 

6 

tt 

tt 

+0.0ii9 

+0.153 

II 

9 

1 

-26.67 

-15.555 

-0.158 

+0 . 329 

0 

1*.8 

-26.9 

Jan 

73 

2 

-23.33 

-ll.llj 

-O.I125 

+0.296 

II 

II 

II 

ooor 

3 

-23.89 

-  9.72 

+0 . 087 

+0.270 

" 

II 

II 

to 

o8oo 

-28.33 

-16.67 

-O.3U3 

+0.287 

II 

II 

5 

-25.555 

-16.11 

-0.11*9 

-0.1*21 

" 

" 

II 

6 

-25.555 

-16.67 

-0.132 

+0.165 

II 

II 

II 

o8oo 

1 

-  2.22 

-Ilt.l4l4 

+0.972 

+0.579 

8.20 

■ 

-16.9 

to 

1700 

2 

-  9.hu 

-13.33 

+0.1885 

+0.1*10 

II 

B 

II 

3 

-12.22 

-11.11 

+1.101* 

+0.605 

" 

" 

II 

+  1.67 

-12.22 

+0.1*76 

+0.362 

II 

B 

II 

5 

-  3.89 

-  8.89 

+0.303 

+0.208 

II 

B 

" 

6 

-  7.22 

-13.89 

+0.212 

+0.51*85 

II 

■ 

II 

1700 

1 

-17.78 

-10.555 

-0.21*0 

+0.3225 

0 

II 

-19.1* 

tc 

2-00 

2 

-15.00 

-10.00 

-0.287 

+0.316 

II 

II 

II 

•3 

-11.67 

-  6.9i* 

+0.01*2 

+0.1*18 

" 

" 

- 

-21.11 

-  9.hh 

-0.326 

+0.133 

II 

B 

II 

C 

y 

-18.33 

-10.555 

-0.196 

-0.376 

II 

B 

" 

c 

-10.555 

-10.28 

-0.075 

+0.077 

II 

B 

II 

2 


Table  XIII.  Snow  depths  as  a  function  of  roof  and  original  snow  covn. 

‘Original  snow  cover’  is  the  snowfall  occurring  in  such  close  proximity  in  time  as  to  be  undissipated  at  the 
time  of  the  first  measurements  of  snow  depth  on  the  roofs  following  its  cessation. 


Original 

Date 

Time 

Snow  depths 

IQMQB 

oi-.o’..l'all 

Snow  depth 

snow  depth 

(cm 

cr.  ■ 

measured 

measured 

‘'l 

*"2 

5 

6.35 

7  Dec  71 

10l*5 

5.08 

3.81 

5.08 

3.81 

1.27 

5.08 

8  Dec  71 

0930 

11 

II 

II 

II 

1.27 

16  Dec  71 

0930 

1.27 

1.27 

1.27 

1.27 

0 

1.27 

II 

1130 

0.635 

0.635 

0.635 

0 

0 

11 

1330 

0 

0 

0 

II 

II 

tl 

7.62 

18  Dec  71 

1115 

0 

1.9 

2.5 

5.715| 

0 

5.715 

to 

to 

5.1* 

6.0 

2.03 

19  Dec  71 

1300 

II 

0.6 

1.9 

1*.13 

It 

3.175 

to 

to 

3.5 

5.7 

6.86 

21  Dec  71 

0800 

5.08 

6.35 

7.62 

8.89 

5.08 

7.62 

M 

11*30 

3.81 

5.08 

II 

7.62 

" 

22  Dec  71 

0800 

2. 51* 

2.5 

3.8 

5.08 

" 

5.08 

to 

to 

3.8 

7.6 

23  Dec  71 

1000 

1.3 

It 

2.5I* 

9 

2.5 

to 

to 

5.0 

5.0 

1300 

II 

11 

11 

11 

H 

2.51* 

2l*  Dec  71 

1010 

It 

11 

11 

II 

It 

1300 

0 

0 

0 

0 

0 

0 

17.78 

31  Dec  71 

1030 

16.51 

20.32 

20.32 

16.51 

9 

17.78 

7.62 

3  Jan  72 

0930 

17.11*5 

19.685 

19.685 

13.97 

5.08 

19.685 

1*  Jan  72 

0800 

15. 2U 

17.78 

15.21* 

8.89 

2. 51* 

16.51 

2.79 

5  Jan  72 

1000 

13.335 

17.11*5 

16.51 

12.7 

5.08 

18.1*15 

10  Jan  72 

1300 

13.0 

17.0 

16.2 

0 

0 

18.1 

1615 

0 

0 

0 

tl 

If 

10.1*1 

21  Jan  72 

1630 

2.51* 

l*.l*l* 

l*.i*l* 

6.35 

6.985 

5.08 

2.51* 

2l*  Jan  72 

1500 

0 

0 

0 

0 

0 

0 

3.56 

28  Jan  72 

1500 

2.51* 

2.5I* 

2.51* 

2.51* 

2. 51* 

2.  51* 

5.81* 

1*  Feb  72 

9 

0.5 

0.5 

3.5 

3.5 

0 

0 

0.76 

7  Feb  72 

It 

0.6 

0.6 

0.6 

0 

- 

11 

to 

to 

to 

2.5  . 

2.5  . 

2.5_ 

1*0.39 

28  Feb  72 

0815 

5.715 

6.35 

10.16 

2.5I* 

2. 51* 

3.81 

29  Feb  72 

t1 

5.6 

6.2 

10.0 

0 

0 

0 

1  Mar  72 

It 

0 

0 

0 

tl 

tl 

It 

6.35 

6  Mar  72 

0800 

2.51* 

3.81 

5.08 

5.08 

3.81 

5.08 

0.25 

7  Mar  72 

It 

It 

2.51* 

2.51* 

2.5I* 

2.51* 

2.51* 

17.78 

15  Mar  72 

0830 

12.7 

12.7 

12.7 

12.7 

trace 

12.7 

151*5 

tt 

It 

13.97 

13.97 

0 

11*.  605 

16  Mar  72 

0800 

11.1*3 

" 

ft 

12.7 

tl 

13.97 

1600 

-  ---1 

tt 

It 

tt 

0 

0 

28 
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Table  XIII  (cont’d). 


Original 

snowfall 

cm 

Date 

anow  depth 

measured 

Time 

snow  depth 
measured 

H 

li 

H 

6.35 

23  Dec  72 
2k  Dec  72 
25  Dec  72 

9 

n 

■ 

■ 

B 

B 

B 

■ 

16.67 

5  Jan  73 

? 

21.59 

23.1*95 

30.1*8 

0 

15.21* 

8  Jan  73 

ti 

" 

II 

II 

15. 21* 

II 

II 

9  Jan  73 

If 

It 

II 

11 

iI*.605 

II 

ii*.605 

10  Jan  73 

II 

II 

II 

II 

13.97 

II 

•1 

11  Jan  73 

ft 

11 

II 

II 

II 

13.97 

12  Jan  73 

II 

11 

II 

12.7 

to 

ll*.0 

II 

13  Jan  73 

It 

" 

II 

M 

13.0 

II 

II 

ll*  Jan  73 

II 

It 

II 

" 

It 

II 

*' 

15  Jan  73 

II 

20.955 

22.86 

29.8U5 

12.7 

II 

12.7 

16  Jan  73 

11 

20.32 

II 

10.16 

II 

II 

17  Jan  73 

It 

" 

II 

27. 9I* 

5.08 

II 

11.1*3 

18  Jan  73 

*' 

5.08 

7.62 

15.21* 

0 

II 

2.51* 

19  Jan  73 

1315 

2.03 

2.51* 

5.08 

II 

0 

22  Jan  73 

•7 

0 

0 

0 

II 

II 

29.1*6 

29  Jan  73 

12l*0 

20.32 

20.32 

20.32 

25.1* 

WBi 

BM 

16  Fet  73 

7 

. .  .  1 

2.51* 

2.5I* 

2. 51* 

2.5I* 

w 

Bm 

the  aluminum  roofs  that  are  the  exception.  The  reason 
for  this  is  that  the  no.  5  (high-slope,  aluminum)  roof 
was  the  roof  with  the  least  snow  cover  and  therefore 
most  able  to  interact  strongly  with  solar  radiation. 
Again,  a  study  of  Table  11  reveals  that  r'  is  least  for 
the  aluminum  roofs  and  therefore,  geometrically  speak¬ 
ing,  these  roofs  had  the  least  ability  to  hold  snow  cover. 

Comparison  of  the  two  types  of  stabUity  rankings  is 
interesting.  While  the  no.  3  (low-slope,  cedar)  roof  had 
the  most  tendency  to  maintain  the  same  temperature 
profile,  when  changes  did  occur,  they  were  more  dras¬ 
tic  changes  than  those  that  occurred  for  the  no.  4 
(high-slope,  asphalt)  roof.  The  implication  is  that  the 
cedar  roofs  changed  profile  less  often  but  more  dras¬ 
tically  than  the  asphalt  roofs  and  that  the  aluminum 
roofs  changed  both  more  often  and  more  drastically 
than  roofs  of  either  asphalt  or  cedar  composition.  This 
result  can  be  understood  by  referring  to  Table  11  and 
noting  that  the  thermal  conductivity  k'  and  thermal 
diffusivity  A '  are  similar  for  asphalt  and  cedar  but  dif¬ 
fer  by  orders  of  magnitude  for  aluminum  compared  to 
either  asphalt  or  cedar. 

The  next  question  of  interest  is  when  did  the  pro¬ 
files  change  and  what  relationships  did  the  changes 
have  to  solar  radiation?  Refer  to  Table  VI.  The  temp¬ 
erature  profiles  of  all  the  roofs  were  most  stable  during 
the  seven  hours  before  dawn  and  became  a  complex 


function  of  the  interaction  between  solar  radiation, 
roof  characteristics,  and  snow  cover  during  the  ten 
hours  from  dawn  to  dusk.  They  then  tended  to  return 
to  their  predawn  status  in  the  seven  hours  after  dusk. 

In  general,  each  roof  tended  to  have  one  characteristic 
profile  in  the  absence  of  solar  radiation  and  another  or 
a  series  of  other  profiles  in  the  presence  of  solar  radia¬ 
tion.  The  ridge,  eaves,  and  low  sections  heated  most 
rapidly  while  the  middle,  under-deck,  and  floor  sections 
responded  most  slowly.  The  no.  5  (high-slope,  alumi¬ 
num)  roof,  as  expected  from  the  results  given  above, 
provided  the  most  extreme  example  of  a  characteristic 
profile  change.  When  solar  radiation  was  applied,  the 
temperature  profile  essentially  reversed  itself  so  that 
the  cooler  sections  in  the  early  morning  became  the 
warmer  sections  at  noon. 

Using  Table  VI  and  ordering  according  to  the  tend¬ 
ency  to  exhibit  this  type  of  profile  change : 
no.  5  (high-slope,  aluminum)  roof 
no.  4  (high-slope,  asphalt)  roof 
no.  3  (low-slope,  cedar)  roof 
no.  2  (low-slope,  aluminum)  roof 
no.  6  (high-slope,  cedar)  roof 
no.  1  (low-slope,  asphalt)  roof 
Except  for  the  no.  3  (low-slope,  cedar)  roof,  this  order 
implies  that  the  controlling  factors  are  k'  and  slope. 

For  each  composition,  the  high  slope  is  more  prone  to 
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Table  XIV.  Snow  coverage  as  a  function  of  roof  and  original  snow  cover. 

“Original  snow  cover"  it  the  snowfall  occurring  in  such  close  proximity  in  time  as  to  be  undissipated  at  the 
time  of  the  first  measurements  of  snow  coverage  on  the  roofs  follo%virig  its  cessation. 


Table  XIV  (cont’d). 


Orifrinal 

snowfall 

(cm) 

Date 

coverage 

measured 

Time 

covera^se 

measured 

H 

COVEFJIG 

’3 

ES  OF  J 

- 

th  POOF 

^5 

’6 

10  Jan  73 

It 

II 

It 

" 

II 

II 

II 

11  Jan  73 

II 

87 

87 

87 

!»0 

20 

II 

12  Jan  73 

If 

If 

It 

II 

II 

13  Jan  73 

II 

II 

If 

II 

II 

lU  Jan  73 

II 

II 

II 

11 

II 

15  Jan  73 

II 

II 

II 

II 

II 

l6  Jan  73 

II 

II 

II 

II 

35 

0 

67 

17  Jan  73 

II 

75 

75 

75 

30 

60 

l8  Jan  73 

67 

67 

67 

0 

17 

19  Jan  73 

1315 

60 

30 

17 

0 

22  Jan  73 

0 

0 

0 

mm 

■i 

29.  h6 

29  Jan  73 

12i*0 

100 

100 

100 

100 

75 

100 

l6  Feb  73 

II 

" 

II 

70 

exhibit  the  characteristic  profile  reversal  (except  for  the 
cedar  roofs).  For  each  slope,  the  compositions  in  order 
of  tendency  to  exhibit  profile  reversal  are  aluminum, 
asphalt,  and  cedar  [except  for  the  no.  3  (low-slope, 
cedar)  roof] .  The  latter  order  is  precisely  the  order 
of  decreasing  values  of  k' . 

Why  does  the  no.  3  (low-slope,  cedar)  roof  appear 
as  an  exception?  Examination  of  Table  VI  reveals  that, 
unlike  the  case  for  the  other  roofs  where  it  is  the  ridge 
and  eaves  sections  that  are  responding,  in  the  case  of 
the  no.  3  (low-slope,  cedar)  roof,  it  is  the  floor  and 
ridge  sections  that  are  responding  most  rapidly.  It  has 
already  been  noted  that  this  roof  retained  too  much 
snow  cover  to  offer  a  true  response  to  solar  radiation. 
However,  it  appears  that  faulty  insulation  in  the  floor 
allowed  heat  from  inside  the  building  to  affect  the  pro¬ 
file  of  the  roof.  As  there  were  only  two  measurements 
of  floor  temperature,  rather  than  one  for  each  roof,  it 
may  be  concluded  that  the  no.  3  (low-slope,  cedar)  roof 
actually  was  the  least  likely  to  exhibit  profile  reversal 
due  to  solar  radiation  and,  just  because  of  its  failure  to 
react  to  solar  radiation,  it  showed  the  weaker  effects  of 
heat  loss  from  the  interior  of  the  building.  In  other 
words,  the  snow  cover  on  top  of  the  roof  insulated  it 
from  solar  radiation  more  effectively  than  the  floor  in¬ 
sulation  from  interior  heat,  which,  like  solar  radiation, 
was  applied  only  during  daylight  hours. 

Common  discussions  in  the  literature’  tend  to  lay 
the  major  blame  for  heating  of  the  middle  section  of  a 
roof  on  warm  attic  air  and  to  discount  the  influence  of 
solar  radiation.  It  is  clear  from  the  results  here  that  the 
true  situation  is  not  so  simple.  Warm  attic  air  and  solar 
radiation  affect  every  roof  and  it  is  the  slope,  roof 
geometry,  and  thermal  conductivity  of  its  covering  that 
determine  which  effect  is  greatest.  In  particular,  for  a 


roof  of  high  slope  covered  with  a  material  of  high  ther¬ 
mal  conductivity,  it  is  solar  radiation  rather  than  warm 
attic  air  that  dictates  the  final  temperature  profile. 
Furthermore,  the  same  is  probably  true  of  any  roof  thai 
does  not  consistently  retain  a  snow  cover  approaching 
30  cm  in  depth. 

What  causes  the  observed  profile  reversal?  The  ridge 
section  responses  are  understandable,  since  the  ridge  is 
the  section  most  exposed  to  solar  radiation  and  warm 
r  ■  :  air.  How  is  it  that  the  eaves  section  heats  more 
rapidly  than  the  middle  section?  This  may  be  evidence 
that  the  eaves  are  heated  from  below  by  warm  air  con¬ 
vection  upwards  on  the  wall  and  collecting  under  the 
eaves.  The  higher  the  slope,  the  greater  the  ability  of 
the  roof  to  retain  warm  air  under  the  eaves  and  the 
greater  the  thermal  conductivity  of  the  roof  composi¬ 
tion,  the  more  rapidly  it  would  heat.  This  collection 
of  warm  air  under  the  eaves  has  been  observed  in  thermo¬ 
grams;  Figure  5  provides  an  example.’ 

Comparative  temperatures  of  roofs  by 
section  -  rankings 

The  temperature  rankings  compare  the  roof  temp¬ 
eratures  for  a  specific  section.  These  data  are  sum¬ 
marized  in  Tables  VII  to  X. 

As  for  the  discussion  of  the  profiles  above,  we  be¬ 
gin  with  the  question  of  simple  stability.  For  which 
section  are  the  roof  temperatures  most  nearly  all  in 
the  same  relationship?  From  Table  VII,  the  answer 
in  order  of  decreasing  stability  is; 
floor 
middle 
under  deck 
ridge 
low 
eaves 
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Figures.  Thermogra^n  oj  eave  section  showing  trapped  heat  under  eave.  (Ref.  7, 
Fig.  9.) 


The  results  for  the  floor  section  are  skewed  somewhat 
because  of  the  fact  that  there  were  only  two  tempera¬ 
ture  measurements  rather  than  one  for  each  roof.  The 
implication  is  then  that  the  middle  section  is  the  sec¬ 
tion  most  characteristic  of  the  roof  as  a  whole. 

Moving  to  Table  Vlll,  the  question  changes  to 
what  section’s  rankings  change  the  least  in  terms  of 
significant  changes?  The  order  becomes: 
floor 
middle 
under  deck 
low 
eaves 
ridge 

This  result  could  have  been  expected,  since  it  was  the 
ridge,  low,  and  eaves  sections  that  were  found  above 
to  respond  most  rapidly  to  the  application  of  solar 
radiation. 

When  do  the  changes  in  temperature  rankings  occur 
for  a  given  section?  The  answer  is  in  Table  IX.  The 
middle  section’s  ranking  remains  nearly  the  same 
throughout  the  day.  At  the  other  extreme,  the  eaves 
section’s  ranking  changes  continually  with  changes  in 
the  meteorological  variables  and  tends  to  reverse  it¬ 
self  ;  the  roofs  that  are  warmest  at  night  tend  to  be¬ 
come  the  coolest  roofs  during  daylight.  In.order  of 
the  tendency  to  show  this  ranking  reversal  with  solar 
cycle : 

eaves 

low 

ridge 

under  deck 

floor 

middle. 

If  the  middle  section  is  the  section  most  character¬ 
istic  of  the  roof  as  a  whole,  its  characteristic  ranking 


is  the  one  that  most  nearly  answers  the  question  of 
which  roof  tends  to  be  coolest  and  which  tends  to  be 
warmest.  The  maximum  ranking  for  the  middle  sec¬ 
tion  is  1421)365.  This  means  that  the  no.  1  (low-slope, 
asphalt),  no.  2  (low-slope,  aluminum),  and  no.  4  (high- 
slope,  asphalt)  roofs  tend  to  respond  alike  and  to  be 
the  coolest  roofs,  while  in  order  of  increasing  tempera¬ 
ture,  the  remaining  roofs  are: 

no.  3  (low-slope,  cedar)  roof 
no.  6  (high-slope,  cedar)  roof 
no.  5  (high-slope,  aluminum)  roof. 

This  implies  that  the  higher  the  slope,  the  greater 
the  solar  heating  of  the  middle  section.  The  fact  that 
the  middle  section  tends  to  be  cooler  for  the  two 
asphalt  roofs  (regardless  of  slope)  implies  that  a  sec¬ 
ondary  effect  due  to  composition  is  involved.  This 
compositional  effect  may  be  due  to  thermal  diffusiv- 
ity  since  this  property  increases  in  the  order  asphalt, 
cedar,  aluminum  like  the  observed  temperature  rank¬ 
ing.  That  is,  the  thermal  diffusivity  of  aluminum  is 
orders  of  magnitude  greater  than  the  thermal  diffus¬ 
ivity  of  the  other  two  materials  and  it  is  the  aluminum 
roofs  that  are  most  different  from  each  other,  and  it 
is  the  asphalt  roofs  that  are  the  most  similar. 

The  effects  of  warm  attic  air  may  be  determined  by 
studying  the  rankings  for  the  under  deck  and  floor 
sections.  Table  VIII  clearly  reveals  that  slope  is  the 
dominant  effect  and  that  the  lower  the  slope,  the 
warmer  the  attic  air.  Subsidiary  compositional  effects 
are  that  the  attic  air  tends  to  be  cooler  under  an  alumi¬ 
num  roof  than  under  a  cedar  roof  but  that  the  slope 
does  not  affect  the  attic  air  temperature  of  an  asphalt 
roof  as  much  as  it  does  roofs  of  the  other  two  mate¬ 
rials.  While  the  no.  1  (low-slope,  asphalt)  roof  has  the 
warmest  attic  air  of  the  low-slope  roofs,  the  no.  4 
(high-slope,  asphalt)  roof  has  the  coolest  attic  air  of 


Table  XV.  Meltwater  volume  collected  subsequent  to  each  snowfall. 

ExproiMd  as  a  volume  and  as  a  fraction  of  the  original  mass  of  snow  cover  on  the  roof  immediately  after 
the  snowfall's  cessation.  Dates  given  are  those  from  cessation  of  the  snowfall  to  dissipation  of  snow  cover 
on  all  roofs  or  the  beginning  of  a  new  snowfall*  whichever  was  applicable. 


Table  XVI.  Degree  of  icing. 

lltlnf  the  uinisiuaiitilative  definition  of  "degree  of  icing.”  viiuel  inspection  wu  used  to  estimste  the  rela¬ 
tive  icing  on  the  test  roofs.  The  visual  data  are  recorded  photographically  in  Figure  4.  Refer  to  Table  I 
for  the  definition  of  1. 


Observation 

date 

Dep 

ree  of 
^2 

icing 

)n  each 

test  rc 

S 

of 

^6 

Fipure  4 
photo  nos. 

•T  Dec  71 

0 

0.75 

0 

0 

0 

0 

15,  17,  19,  22 

20  Dec  71 

0 

0 

0 

0 

0 

0.25 

31,  32,  31,  36 

21  Dec  71 

0 

1.00 

0 

0 

0 

0 

39,  lO,  ll,  12 

23  Dec  71 

0.50 

0.50 

0.50 

0.25 

0.25 

0.25 

17,  19 

3  Jan  72 

0 

0.50 

0 

0 

0 

0 

50,  51,  52,  53 

Jan  72 

0 

1.00 

0.75 

0.25 

0.50 

0.25 

6^,  65,  66,  6T 

5  Jan  72 

1.00 

1.00 

1.00 

0 

0.50 

0.50 

69,  70,  71,  72 

7  Jan  72 

0 

1.00 

0.25 

0 

0.25 

0 

75,  ■’7,  78,  •’9 

10  Jan  72 

1.00 

1.00 

1.00 

0 

0 

0 

81,  82,  81 

28  Feb  72 

0.50 

0.50 

0.50 

0.25 

0.50 

0.50 

96,  98,  99,  102, 
101,  105,  10” 

l6  Mar  72 

0 

0 

0 

0 

0.25 

0.25 

118,  119,  120,121 

H 

0.75 

0.75 

0.75 

0 

0 

0 

127,  128,  129,130 

the  high-slope  roofs  and,  in  fact,  is  more  like  its  low- 
slope  counterpart  than  the  other  high-slope  roofs. 

Since  the  middle  section’s  maximum  ranking  is  not 
at  all  similar  to  the  under-deck  section’s  maximum 
ranking,  1421)365  as  compared  with  564123.  respec¬ 
tively,  it  is  clear  that  solar  radiation  rather  than  warm 
attic  mr  had  the  most  effect  on  the  roofs  as  a  whole. 

What  effect  is  dominant  in  producing  the  ranking 
reversals?  Compositional  effects  must  be  separated 
from  slope  effects.  Refer  to  Tables  X  and  XI.  In 
Table  X,  all  of  the  possible  equal  temperature  cases 
are  listed.  Study  of  the  percentage  of  the  total  test 
hours  for  which  each  case  was  observed  for  each  sec¬ 
tion  (exclusive  of  the  floor  section)  reveals; 

1)  Only  the  middle  and  low  sections  exhibited  any 
of  the  possible  equal  temperature  cases  for  more  than 
50%  of  the  test  hours. 

2)  The  equalities  (41)  and  (21)  were  the  only  equal¬ 
ities  observed  for  more  than  50%  of  the  test  hours. 

3)  For  all  sections,  either  (41)  or  (2 1)  was  the 
equality  most  often  observed. 

These  results  agree  with  previous  indications  that,  for 
asphalt,  compositional  effects  are  dominant  and,  for 
aluminum  and  cedar,  slope  and  geometrical  effects 
are  dominant.  Table  XI  shows  the  data  analyzed  from 
the  opposite  point  of  view  and  shows  the  percentage 
of  test  hours  for  which  a  deflnitive  slope  dependence 
was  observed.  This  analysis  implies: 

1)  Attic  air  temperature  increases  with  slope. 

2)  Ridge  temperature  decreases  with  slope. 

3)  Middle  and  low  sections  increase  in  temperature 


with  slope  but  do  not  respond  to  a  change  in  slope 
as  readily  as  the  ridge,  under-deck,  and  floor  sections. 

4)  The  eaves  section  shows  a  stronger  depiendence 
upon  composition  than  on  slope. 

Heat  flow  and  accumulation 

For  the  two  test  days  that  were  compared,  one  had 
much  less  solar  radiation  than  the  other.  In  addition, 
for  the  test  day  with  the  most  solar  radiation,  the 
snow  cover  was  deepest  but  least  uniform,  with  the 
no.  5  (high-slope,  aluminum)  roof  being  bare  and  the 
eaves  of  the  other  roofs  exposed.  For  the  day  of  least 
solar  radiation,  the  snow  covet  was  thin  but  uniform 
and  all  sections  of  every  roof  were  covered.  (Refer 
to  Table  Xll.) 

For  the  day  with  low  solar  radiation,  the  heat  bal¬ 
ance  of  all  the  roofs  for  both  the  middle  and  eave 
sections  did  not  change  signiflcantly  throughout  the 
solar  cycle.  The  no.  1  (low-slope,  asphalt)  and  no.  5 
(high-slope,  aluminum)  roofs  tended  to  lose  heat  and 
the  no.  2  (low-slope,  aluminum),  no.  4  (high-slope 
asphalt),  and  no.  6  (high-slope,  cedar)  roofs  tended  to 
gain  heat  throughout  the  day  for  the  eaves  section. 

For  the  middle  section,  all  roofs  except  no.  1  (low- 
slope,  asphalt)  tended  to  gain  heat  during  the  daylight. 

For  the  day  with  high  solar  radiation,  the  heat  bal¬ 
ance  changed  markedly  for  the  eave  section  (exhibit¬ 
ing  an  increase  in  heat  accumulation  for  all  roofs)  on 
all  roofs  but  did  not  change  much  for  the  middle  sec¬ 
tion  except  in  the  case  of  the  no.  5  (high-slope,  alumi¬ 
num)  roof.  For  the  no.  5  (high-slope,  aluminum)  roof. 
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Figure  6a. 

Figure  6.  Change  in  coverage  of  nxifs  with  time. 
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both  the  eave  and  the  middle  sections  exhibited  a 
marked  rise  in  heat  accumulation  during  daylight.  This 
was  probably  due  to  the  lack  of  snow  cover  of  this  roof 
for  this  test  day. 

In  general,  these  results  agree  with  the  results  ob¬ 
tained  from  the  thermocouple  data. 

Snow  depth,  coverage,  and  meltwater 

The  snow  depth  was  generally  least  for  the  high- 
slope  roofs  and,  for  a  given  slope,  increased  with  r', 
the  size  of  the  butt  end  of  the  shingles.  The  exception 
to  this  was  the  no.  2  (low-slope,  aluminum)  roof, 
which  had  the  greatest  tendency  to  form  ice  dams  at 
the  eave  and  thereby  hold  a  deeper  snow  cover.  Refer 
to  Table  Xlll. 

For  a  given  roof,  the  eave  section  tended  to  become 
bare  first  followed  by  the  ridge  section.  The  no.  5 
(high-slope,  aluminum)  roof  often  cleared  days  sooner 
than  the  other  roofs.  This  can  be  seen  in  the  sequence 
of  phonographs  shown  in  Figure  6  and  in  the  coverage 
data  detailed  in  Table  XIV. 

The  meltwater  as  a  fraction  of  original  snow  cover 
averaged  for  the  18  snowfalls  recorded  in  Table  XV 
decreased  in  the  order: 

no.  5  (high-slope,  aluminum)  roof 
no.  4  (high-slope,  asphalt)  roof 
no.  2  (low-slope,  aluminum)  roof 
no.  1  (low-slope,  asphalt)  roof 
no.  3  (low-slope,  cedar)  roof 
no.  6  (high-slope,  cedar)  roof. 


In  comparison,  the  total  meltwater  volume  averaged 
for  the  same  snowfalls  decreased  in  the  same  order  ex¬ 
cept  that  the  no.  5  (high-slope,  aluminum)  and  no.  6 
(high-slope,  cedar)  roofs  exchanged  positions.  While 
the  no.  5  (high-slope,  aluminum)  roof  produced  the 
least  meltwater  volume  pet  snowfall,  a  greater  fraction 
of  its  snow  cover  ran  off  as  meltwater  following  the 
snowfall’s  cessation.  The  no.  6  (high-slope,  cedar)  roof 
produced  the  most  meltwater  volume  but  the  least 
fraction  of  its  snow  cover  ran  off  as  meltwater. 

If  the  results  above  are  compared  with  the  results 
for  the  temperature  ranking  for  the  middle  section,  the 
data  seem  to  say  that  the  cooler  the  roof  the  more  melt¬ 
water  it  produces.  This  apparent  contradiction  is 
erased  when  the  effects  of  slope  are  considered.  The 
meltwater  produced  by  the  high-slope  roofs  tends  to 
run  off  more  rapidly  than  that  produced  by  the  low- 
slope  roofs,  which  freeze  more  of  their  meltwater  into 
ice  dams  because  of  the  slower  rate  of  flow. 

Degree  of  icing 

Table  XVI  records  the  observed  degree  of  icing,  /. 

For  the  12  observation  times,  the  average  value  of  /  de¬ 
creased  in  the  order; 

no.  2  (low-slope,  aluminum)  roof 
no.  3  (low-slope,  cedar)  roof 
no.  1  (low-slope,  asphalt)  roof 
no.  5  (high-slope,  aluminum)  roof 
no.  6  (high-slope,  cedar)  roof 
no.  4  (high-slope,  asphalt)  roof. 
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This  result  implies  that  the  higher  the  slope,  the  less 
likely  the  roof  was  to  ice;  and  for  a  given  slope,  alumi¬ 
num  was  the  most  likely  and  asphalt  the  least  likely  to 
ice. 

Discussion  of  the  icing  problem 

If  the  semi-quantitative  indicator  /  could  be  expressed 
as  a  quantitative  function  of  the  relevant  variables,  the 
following  assumptions  might  reasonably  be  stated: 


D/^CT’M-rE)  if  7’m>7’e 

(3) 

2)  /«  1/K 

(4) 

3)  lech 

(5) 

4)  lay 

(6) 

5)  la  1/a 

(7) 

6)  lar' 

(8) 

7)  Icck' 

(9) 

8)  laA' 

(10) 

9)  las'. 

(11) 

All  of  these  assumptions  are  supported  to  some  extent 
by  the  results  discussed  previously. 

If  it  is  assumed  that  each  of  the  nine  assumptions  is 
a  valid  statement  about  icing,  which  of  the  nine  effects 
is  most  important?  The  answer  is  clearly  assumption 
(S).  The  slope  of  the  roof  dominates  its  icing  behavior. 

It  does  this  by  controlling  h,  y,  and  V.  If  all  the  snow 
leaves  the  roof  as  meltwater  or  faUen  snow,  none  will 
be  left  to  freeze  into  ice.  The  assumptions  concerning 
the  effects  of  k'  A'  5' and  r '  are  the  next  most  im¬ 
portant  statements,  r'  determines  the  ability  of  the 
roof  to  hold  snow  cover  for  a  given  slope;  i.e.,  a  stepped 
slope  can  hold  snow  more  readily  than  a  smooth  one  and 
the  greater  the  size  of  the  steps,  the  greater  the  ability  to 
hold  snow,  k' ,A\ and  S'  determine  the  thermal  inter¬ 
action  of  the  roof  covering  with  the  snow  cover,  solar 
radiation,  heat  leaking  through  the  attic  floor,  and  heat 
trapped  under  the  eaves.  They  determine  the  value  of 
^M'^E  profile  of  the  roof. 

The  nine  assumptions  above  could  therefore  be  sum¬ 
marized  by 

/orm'A'/a  (12) 

and 


[T„ik'X^').T^ik'^A')]hilla.T')y{llay) 

“  y(a,k',S'^') 

(13) 

where  eq  12  indicates  the  dependence  of  icing  on  the 
roof  characteristics  directly  and  eq  13  depicts  the  man¬ 
ner  in  which  the  roof  characteristics  act  to  control  the 
icing.  Equations  12  and  13  do  not  represent  quantita¬ 
tive  results;  they  represent  the  simplest  form  that  quan¬ 
titative  results  could  take.  It  is  therefore  clear  that 
each  of  the  roof  characteristics  acts  in  an  interrelated 
manner  to  affect  V,  h,  y,  and  the  temperature  profile 
of  the  roof,  the  conditions  that  lead  to  the  formation 
of  ice.  While  the  slope  is  the  most  important  character¬ 
istic,  the  other  roof  characteristics  may,  if  they  assume 
extreme  values,  enhance  or  counter  the  effect  of  the 
slope.  Furthermore,  it  is  clear  that,  as  the  size  and 
frequency  of  snowfalls  increase,  the  roles  of  the  ther¬ 
mal  characteristics  will  become  more  important  be¬ 
cause,  once  any  ice  has  formed,  it  acts  as  a  barrier  to 
the  flow  of  falling  snow  and  meltwater,  thereby  effec¬ 
tively  reducing  the  ability  of  the  slope  to  prevent  icing. 
Thus,  for  flat  or  nearly  flat  roofs  (such  as  the  pro¬ 
tected  membrane  rooO  the  importance  of  the  insula¬ 
tion  and  thermal  and  optical  properties  of  the  surface 
covering  cannot  be  over-emphasized.  The  need  for 
increased  insulation  in  protected  membrane  roofs  has 
been  emphasized  elsewhere.'  The  results  here  indi¬ 
cate  that  the  choice  of  surface  covering  must  receive 
equal  attention. 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 
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